Published by Institute of Physics Publishing for SISSA/ISAS 

Received: January 29, 2009 
Accepted: January 29, 2009 



OS 

o 
o 



(N 



The origin of 'Great Walls' 



rgei F. Shandarin 

Department of Physics and Astronomy, UnroersUy of Kansas, KS 66045, U.S.A. 



OS ■ E-mail: 



sergeiOku . edu 

o 



43 ■ Abstract: A new semi-analytical model that explains the formation and sizes of the 'great 
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o 

H 



O 
L\3 



the model is the Zel'dovich approximation it has been used in a new way very different 
from the previous studies. Instead of traditional approach that evaluates the nonlinear 
density field it has been utilized for identification of the regions in Lagrangian space that 
after the mapping to real or redshift space (depending on the kind of structure is studied) 
end up in the regions where shell-crossing occurs. The set of these regions in Lagrangian 
space form the progenitor of the structure and after the mapping it determines the pattern 
, of the structure in real or redshift space. The particle trajectories have crossed in such 

"n^ \ regions and the mapping is no longer unique there. The progenitor after mapping makes 

only one stream in the multi-stream flow regions therefore it does not comprise all the 
mass. Nevertheless, it approximately retains the shape of the structure. The progenitor 
of the structure in real space is determined by the linear density field along with two non- 
Gaussian fields derived from the initial potential. Its shape in Eulerian space is also affected 
^ ■ by the displacement field. The progenitor of the structure in redshift space also depends on 

these fields but in addition it is strongly affected by two anisotropic fields that determine 
the pattern of great walls as well as their huge sizes. All the fields used in the mappings 
are derived from the linear potential smoothed at the current scale of nonlinearity which 
is R n i = 2.7 h Mpc for the adopted parameters of the ACDM universe normalized to 
erg = 0.8. The model predicts the existence of walls with sizes significantly greater than 
500 h -1 Mpc that may be found in sufficiently large redshift surveys. 
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1. Introduction 

The first indications of the existence of superclusters of galaxies bridging the clusters of 
galaxies [l], ^] was confirmed by the CfA redshift survey ||. In the following two decades 
the discovery of huge concentrations of galaxies spanning over 200 h _1 Mpc were reported 
(e.g., 0, |5|). They were dubbed "great walls" due to their sizes and geometry in redshift 
space. The current record belongs to "a Sloan Great Wall of galaxies 1.37 billion light 
years long, 80% longer than the Great Wall discovered by Geller and Huchra and therefore 
the largest observed structure in the universe" ||. Taking into account that the SDSS 
redshift survey is much larger and deeper than CfA one may wonder whether the Sloan 
Great Wall will remain the greatest wall or even greater walls will be found in the future 
redshift surveys and if it is so how large it may be. The answer to this question may 
depend to certain extent on the exact definition of the walls. For instance, some definitions 
may result in the percolating system of filaments and walls that would span throughout 
the whole volume. However even in this case the essential constituents can be probably 
identified and measured. We show that by restricting the analysis to the issues of overall 
geometry and scale of the largest walls one can alleviate some of the problems of this kind 
and make some progress. In particular, the details of the density distribution within the 
walls, such as the accurate positions and masses of halos become less important if the scales 
of interest exceeds ~10 h Mpc. 

We begin with a brief discussion of the structure in real space. The scale separating 
the nonlinear regime of the gravitational growth from linear or mildly nonlinear regimes 



is roughly around 5 h Mpc which is close to the galaxy correlation scale. It seems to be 
by far too small to be directly relevant to the scale of filaments some of them are a least 
70 - 100 h — 1 Mpc long. Bond, Kofman and Pogosyan [|| (hereafter BKP) suggested that 
the filamentary network in real space with a typical scale of ~30 h _1 Mpc can be explained 
by invoking the correlation bridges between relatively high and therefore rare peaks in the 
linear density field filtered on scales greater than such that a p (Rb) < 1. Unfortunately, 
the BKP model has not provided a framework for a quantitative evaluation of the density 
contrast in the filaments bridging the clusters of galaxies from the enhancement of the 
conditional correlation function between two or more peaks. Besides, the gravitational 
growth of the structure in the universe is a deterministic process if the linear perturbation 
field is specified. Therefore, one may seek a deterministic relation between the initial field 
and the final structure at least in cosmological simulations where the full information is 
available. In addition, the cluster analysis of the nonlinear dark matter density field in 
real space obtained by N-body simulation in the ACDM model [jjj revealed the filaments 
with lengths up to ~100 h _1 Mpc p§ i.e., significantly greater than 30 h 1 Mpc. A similar 
analysis of the mock galaxy catalogs || demonstrated the presence of even longer filaments 
in redshift space, up to ~ 150 h Mpc [1C]. Both numbers are almost certainly affected 
by the simulation box sizes (240 h _1 Mpc in the former and 346 h _1 Mpc in the latter case) 
and could be even longer in larger boxes. The simulated dark matter density fields and 
especially mock galaxy catalogs look similar to the observed distribution of galaxies as well 
as they have similar statistics of various kinds that makes the model viable. However, the 
emergence of the nonlinear structures (great walls) ~ 50 - 100 times greater than the scale 
of nonlinearity remans unexplained even in the cosmological N-body simulations where the 
full information about the both initial conditions and formed structure is available. 

It has been realized for a long time that the structures in real space and observed in 
redshift space have different pdfs, power spectra, correlation functions, and higher order 
moments (e.g., [11, [l3|, [14|, [15L |H]]). The redshift-space correlation function has an 



oblate shape on scales greater than a few h _1 Mpc indicating the presence of anisotropic 
structures in redshift space flattened in the radial direction. However, the measurements 
do not go beyond ~ 30 h Mpc where the correlation function drops to ~ 0.1 and the 
signal drowns in noise fll~6|| . Theory has not made any specific predictions concerning the 
structure in the universe on scales greater than 100 h Mpc therefore the conventional 
wisdom is that the perturbations on such scales must be in the linear regime. 

The 'finger of God' is the most conspicuous and best understood effect in redshift 
space. However, it is a relatively small scale effect and cannot explain the coherence of the 
great walls over a few hundred h _1 Mpc. It has been also demonstrated in [17, 18] that the 
size of the structures in the redshift space is considerably greater than that of the parent 
structures in real space. In addition, the authors emphasized a characteristic circular 
pattern in redshift space and the increased spacing between the structures in the redshift 
direction. However, they neither mentioned the enlargement of sizes of the structures in 
the transverse direction nor provided explanation to their observations. 

The presence of walls up to 150 h —1 Mpc in length in redshift space in three-dimensional 
simulations based on the adhesion approximation was emphasized in IS]. The authors 



pointed out that the correlation function was weak or even negative at these scales. They 
also observed that ~ 100 h _1 Mpc walls grow from several favorably aligned high-density 
peaks in the linear density field, each coherent over ~ 20 h _1 Mpc. This seems to be 
a factual description of the process, however similarly to the BKP model it requires a 
dynamical model different from the linear theory of gravitational instability. The major 
problem with the structures spanning over ~ 300 h _1 Mpc arises from a lack of a natural 
scale of this magnitude in the linear density field. If the field is smoothed with say 100 
h _1 Mpc or greater scale than its amplitude becomes so low that such a field would remain 
almost perfectly linear at the present time. 

We view the formation of the structure as a continuous mapping from the initial prac- 
tically uniform state in Lagrangian space to the final highly inhomogeneous and anisotropic 
distribution of mass on scales up to ~ 100 h _1 Mpc in Eulerian space followed by another 
mapping to redshift space where the largest structures become greater than 300 h _1 Mpc. 
These mappings involve the transport of mass (although not physical in the case of the 
mapping to redshift space), however the rms displacement of mass elements is only about 
15 h _1 Mpc and therefore it cannot produce the structures of needed sizes by itself. Thus, 
the certain initial fields i.e., the fields present in the linear stage and playing an important 
role in the nonlinear dynamics must have scales in the range of 100 h _1 Mpc or even greater. 
The major goal of this paper is to show that the fields having needed properties exist in 
Lagrangian space. The outcome of their combined effect is the definition of two parent 
structures that are the progenitors of the large-scale structures: one in real and the other 
in redshift space. Being closely related two progenitors are still quite distinct. The prop- 
erties of the progenitors are determined by a combination of Gaussian and non-Gaussian 
fields derived from the initial gravitational potential smoothed with a Gaussian filter at the 
scale of nonlinearity R n \ defined by the equation < 5f in (R n \) >= 1, where 5 = 5p/p. In the 
standard ACDM model normalized to erg ~ 0.8 this scale is R Q \ = 2.7 h _1 Mpc. The choice 
of the filter and filtering scale for 2D illustrations approximately corresponds to an optimal 
choice of the filter and scale in 3D simulations that studied the accuracy of the Truncated 
Zel'dovich Approximation (TZA) in a set of power law models in the Einstinen-de Sitter 
model [22 1 . Some adjustment will be probably needed for the ACDM model in 3D. Based 
on that study I would conservatively guess that the change in R Q \ would not be greater 
than 30%. No fields smoothed at greater scales are involved in the mapping. The structure 
resulting from the mapping of the progenitors cannot and is not supposed to reproduce 
accurately the mass or galaxy number densities, however the overall shape and sizes of the 
structure in real space and great walls in redshift space must be depicted quite well. 

We use TZA as a mapping device for both mapping [20, [2^, ^l], |15[| . In this paper 
we will concentrate on the lengths of the progenitors in the Lagrangian space and their 
transformation caused by the mapping to Eulerian and redshift space. 

We illustrate the model by two-dimensional plots. Although no two-dimensional model 
can reproduce all the properties of a three-dimensional system it may provide a reasonably 
good description of major ideas. In order to make the appearance of illustrations more 
realistic we use the linear power spectrum P 2D (k) = kP 3D (k), where P 3D (k) is the power 
spectrum of the linear density perturbations in the ACDM model with parameters h = 0.7, 



n m = 0.3, n b = 0.047, n = 1. We use the BBKS transfer function @ with the T- 
parameter given by This choice of two-dimensional power spectrum retains all the 

moments < k 3 > of the three-dimensional spectrum and therefore all the scales of the 
two-dimensional field determined by the ratio of two moments (e.g., the scale of peaks 
R 2 ~< k 2 > / < fc 4 >) remain similar to that of the three-dimensional field except a small 
factor \/3/2 = 1.2. The two-dimensional fields were generated in 1024 h _1 Mpc box on 
1024 2 mesh. 

The rest of the paper is organized as follows. Section 2 outlines the dynamical model 
and fields involved in both mappings. Section 3 provides two dimensional illustrations that 
elucidate the main concepts of the model, the axes in the figures showing various fields 
are given in h _1 Mpc. Section 4 discuss the differences between 2D and 3D cases, Sec. 5 
summarizes the results and outline the further problems, finally the appendix gives the 
joint density probability function of three invariants of the deformation tensor that are 
used in the mappings. 

2. Model 

The dominant feature of the process of the structure formation is usually dubbed as the 
hierarchical clustering. It simply means that small halos generally form earlier than the 
larger ones. Thus, more massive halos are formed by consecutive merging of smaller halos. 
The merging process effectively stops on the scale of R n \ or on a little greater scale. It 
is a well established fact that he formation of the structures on larger scales is mainly 



determined by the long wave part of the initial spectrum with k < [25]. [25] also 
pointed out that the changes in the high-frequency (k > K~^) components of the initial 
conditions did not affect much the evolution of long waves. A direct comparison of the 
nonlinear density fields obtained in the N-body simulations of models with and without 
high-frequency components (Truncated Zel'dovich Approximation) in the initial conditions 
has also shown a very good agreement 22, 21 1. In general, the Zel'dovich approximation 
"gives a reasonable description of the overall pattern of N-body simulations" Q. 

Using the comoving coordinates x = v/a{t) where a(t) is the scale factor describing the 
uniform expansion of the universe the Zel'dovich approximation p6[ (see also pTj) becomes 



x i (c i ,t)=q i +D(t)s i (q), (2.1) 

where vector q gives the Lagrangian (initial) position of a fluid element, x is the Eulireian 
position of the fluid element at time t, D(t) describes the linear growing mode of fluc- 
tuations, Sj(q) = —d<&/dqi is the displacement field and <3?(q) is the linear gravitational 
potential field. Thus, our dynamical model is the Zel'dovich approximation applied to the 
initial fluctuation field smoothed over the nonlinear scale R D \ which equals 2.7 h _1 Mpc in 
the ACDM model normalized to as = 0.8. 

Instead of traditional approach that tries to evaluate the evolution of density field we 
suggest to consider its inverse i.e., the evolution of specific volume V = p~ x . The specific 
volume is more suitable quantity for our purpose because it allows an easy separation of 
the effects related to various fields. We would like to remind that we do not use this model 



for accurate predicting of the nonlinear density field but only for approximate identification 
of high density regions. The specific volume is given by the Jacobian 



<V> 



(2.2) 



The growth of density corresponds to the decrease of the specific volume. Some fluid 
elements may collapse to zero volume and at later time their volumes predicted by eq. 
|2.2| become even negative however this is an easy problem. Negativeness of the density 
or specific volume simply reflects the fact that such a fluid element turns inside out and 
the positiveness of the density can be easily restored by taking the absolute value of the 
Jacobian. It also indicates that the trajectories of fluid elements are crossed and the 
mapping in no longer unique. This creates a much more serious problem that does not 
have a simple solution even in the frame of the Zel'dovich approximation (apart from 
breaking the approximation within the regions where shell-crossing occurs) mainly because 
the problem becomes nonlocal. The set of fluid elements that collapse by the present time 
form a structure in Lagrangian space, we shall call it the progenitor or parent structure. 



The parent structure mapped to Eulerian space by means of eq. |2,l| forms the pattern of 
the large-scale structure. It is worth stressing that the parent elements do not represent 
the total mass contained in the large-scale structure however they approximately mark 
the regions of highest density. The situation becomes more complex if we recall that the 
large-scale structure forms from highly nonuniform state where practically all the mass 
is already in the form of virialized halos. However, as the previous studies showed ( 
e.g.,p7], 25, pOl, 22, ETJ) the overall distribution of the halo number density corresponds to 



this model quite well. 

The Jacobian can be expressed in terms of the invariants (Ii, 12,13) ° r eigen values 
(Ai, A2, A3) of the deformation tensor dij = —dsi/dqj = d 2 <& / dqidqj 



■4(q>*) 



1 - D(t)h(q) + D 2 (t)I 2 (cD - £> 3 (t)/ 3 (q) 
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(2.3) 



where 



^11 + ^22 + ^33, h = M u + M 22 + M33, h 



dn d\2 dn 

d\2 d22 d>23 
dl3 d-23 



(2.4) 



and My are the minors of the corresponding elements of the determinant I3. The eigen 
values are the solutions of the characteristic equation 



A 



JiA 2 + / 2 A-/ 3 



0. 



(2.5) 



It always has three real solutions because the mapping is generated by the potential vector 
field Sj(q) and therefore the deformation tensor is symmetric dij = dji = d 2 Q /dqidqj] it 
is convenient to assume the solutions to be ordered at every point Ai(q) > A 2 (q) and 
A 2 (q) > A 3 (q). 
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Equation J x (q, t) = has obviously three solutions for D as well: Di = 1/Aj. If any 
of them lies in the interval < -Dj(q) < 1 we consider the fluid element with Lagrangian 
coordinates q to belong the progenitor of the structure. Condition < Ai(q) < 1 is 
equivalent to the above condition for -Dj(q) but it is less general and cannot be used for 
the mapping to redshift space where it is not of potential type. 



In order to describe the mapping to redshift space (Q) we supplement eq. 2A with 
additional term (see e.g., JIB]) 



Ci = x i + /£>s fc (q)n fc (x)ni(x) = & + Dsi(q) + fDs k (q)n k (x)ni(x.) (2.6) 
where x is given by eq. 2J., / = ci log -D/d log a « fi ? ° n 6 for the present time |y]] or a better 



approximation given in |p8[ . The unit vector rii = xijx is directed along the line of sight. 
We assume summation over repeated indices throughout the paper. The peculiar velocity 
= afDsi and therefore the displacement of a galaxy from comoving Eulerian position, 
Xi to redshift space position, is v^rik/{aH) = jDsf-n^ along the line of sight. The 
distances in Lagrangian, comoving Eulerian and redshift space are measured in h _1 Mpc. 




Figure 1: Lagrangian space. Regions S(q) > 1 or equvalently 1 — /i(q) < are shown in green. 
The specific volume is equal to the corresponding Jacobian 

— (||). 

The evaluation of in terms of the initial fields is a quite laborious task. Fortunately it 
allows a considerable simplification in the limit of |L*s| <C q. Physically it means that the 
results are not reliable within 50 h _1 Mpc or so from the center but it is not very important 



since we are interested in much greater scales. The components of the Jacobian become 
(to the zeroth order in \Ds\/q) 

« Sij - D(dij + fhih k dj k ), (2.8) 

where n = q/g. Mapping to redshift space is neither irrotational nor solenodial which 
means that Jacobian has no symmetry. Similarly to the Jacobian of the mapping to 
Eulerian space the determinant can be expressed in terms of the invariants of the initial 
fields and the radial components of two tensors d qq and M qq 

J c = l-Dh + D\l + f)I 2 - D\l + f)h ~ fDd qq - fD 2 M qq , (2.9) 

where d qq = hihjdij and M qq = hihjMij are the radial components of the deformation 
tensor and tensor made of the minors of the determinant ^3 . The invariants are statistically 
isotropic fields although two of them (/2(q) and -^(q)) are non-Gaussian. The last two 
terms depend on highly anisotropic fields: one d qq (q) is Gaussian and the other M qq (q) is 
non-Gaussian. 



3. Two-dimensional illustration 
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Figure 2: Lagrangian space. Red regions 1 — 1\ (q) + 12 (q) < 0, green contours S(q) = 1 correspond 
to the boundaries of the regions in fig. [|. 

Major although not all features of the model can be illustrated by a two-dimensional 
model. Although the main goal of the paper is a study of the mapping to redshift space we 



briefly discuss the mapping to real space and show the importance of non-Gaussian fields 
for building the large-scale structure. We shall use a convenient normalization of function 
D{t) such that at the present time D{to) = 1 and the variance of the linear density contrast 
field erf = 1, parameter / ~ 0.8 for the chosen cosmological model. 

3.1 Mapping to real space 

The Jacobian of the mapping to Eulerian space in two dimensions depends on two statis- 
tically isotropic and homogeneous fields: one Gaussian Ii(q) and one non-Gaussian ^(q) 

= J x = 1 - Dh + D 2 I 2 . (3.1) 

The initial density contrast field Si n = 5p/p = D{ti n )I\ is assumed to be Gaussian. Let 
us begin with the linear approximation V x / <V > = 1 — DI±. Linear theory formally 
predicts that the peaks with S > 1 shrink to the regions having formally negative vol- 
umes (and therefore negative density) by the present time. Figure [l] shows an exam- 
ple of the progenitor of the structure in Lagrangian space predicted by the linear theory 
V x (to, q)l < V > = 1 — D(to)Ii(q) = 1 — I\(q) < 0. Although the initial power spectrum is 
filtered with Gaussian windows exp(— k 2 R 2 /2) with Rf = 2.7 h _1 Mpc some peaks S > 1 
have quite large lengths reaching ~ 30 h Mpc. 

The full nonlinear models (t , ?) = I- D{to)I l (q)+D 2 {to)I 2 {q) = I- h(q) + I 2 (q) < 
predicts the progenitor of the structure shown in Figure |2| in red. The green contours 
show the boundaries of the linear progenitor shown in fig.||, thus one can see both the 
similarities and differences of two fields. In order to make an eye ball estimate of sizes 
easier we plot the grid on the background. The similarity of two fields is obvious but it 
is not unexpected because both fields depend on I\(q) and therefore strongly correlated. 
The effect of non-Gaussian 12-field is quite remarkable. The red non-Gaussian regions are 
considerably longer than the linear density peaks. They also typically comprise several 
green regions therefore non-Gaussian field I 2 provides "bridges" linking the peaks of linear 
density contrast field. This is in a good qualitative agreement with the observations made 
by Weinberg and Gunn (ll| who emphasized that the large structures usually consist of 
several peaks in the linear density field. There is also no contradiction to the role of the 
correlation between linear density peaks in the BKP model H, although the the higher 
peaks of the field smoothed with greater scale is assumed in the BKP model. The green 
contours totally comprise 15.9% of the area. It is remarkable that the red regions where 
the nonlinear condition 1 — I%(q) + h(q) < is satisfied occupy slightly less area (15.3% ) 
than the green contours. Despite the greater lengths of the nonlinear progenitor reaching 
over 100 h _1 Mpc red regions are generally slimmer than their counterparts in the linear 
progenitor. The resolution of fig.|2| is far too low to reveal this completely. However, it is 
obvious in fig.|| where the central part of fig.§ is zoomed in. In addition, small isolated 
regions of the linear progenitor disappear from the nonlinear one (empty green contours). 

Figure § shows the structure in Eulerian space (red) along with the progenitor in 
Lagrangian space (green). The mapping makes the structure considerably slimmer than 
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Figure 3: Lagrangian space, the central 200 h 1 Mpc box of fig.||. Red regions 1 — Ji(q)+/2(q) < 0, 
green contours 6(q) = 1 correspond to the boundaries of the regions in fig]]} 



its progenitor in Lagrangian space although it is probably thicker than in reality 
Two points are worth stressing. First, neither lengths of the filaments no topology visibly 
change due to mapping itself. And second, the filaments do not move much, however those 
that moved were are displaced mostly in the transverse direction while the longitudinal 
displacement is significantly smaller. It is also worth recalling that the r.m.s. displacement 
of fluid elements from the Lagrangian to Eulerian positions is about 15 h _1 Mpc which 
greater than the thickness of the filaments but considerably smaller than their lengths. 
The mapping to real space also involves statistically isotropic Gaussian vector field Sj(q) = 
-d$/d qi . 

3.2 Mapping to redshift space 



We obtain the Jacobian of the mapping to redshift space in two dimensions from eq. |2Jj 



J c (q,t) = l-Dh+ D 2 (l + f)I 2 - fDd, 



qq, 



(3.2) 



where I\ = dn + d 2 2 and I 2 = dud 2 2 



d 2 



12 

^2 



are the invariants of the deformation tensor in 



two dimensions, d qq = hfdn + 2nin2<ii2 + ^2^22 is the radial component of the deformation 
tensor. Compared to the mapping to Eulerian space the term proportional to I2 is boosted 
by factor 1 + /. The Jacobian contains a highly anisotropic field d qq (q) in the sense that 
the contours are significantly longer in the transverse direction than in radial direction 
as a result the field acquires a characteristic circular pattern (fig.[|). The cause of this 
anisotropy becomes much more obvious if one considers the small angle approximation 
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Figure 4: Eulerian space. The large-scale structure in Eulerian space (red) and the progenitor of 
the structure in Lagrangian space (green). 



12]. Now, let us consider the direction corresponding to gi-axis. In this case n = (1,0) 
therefore d qq = d\\{q) = d 2 Q(q)/dqf. Differentiating potential two times with respect to 
qi multiplies the power spectrum by kf (P^ xi (k) = kfP^(k)) that effectively reduces the 
characteristic scale of the resulting field along q\ direction while the scale along q<i changes 
much less. In the case of arbitrary direction the reduction of the scale along radial direction 
is much stronger than in the transverse direction that generates a distinct circular pattern. 

The mapping to redshift space also involves the radial component of the velocity field 
v q( c l) = fDs q (q) = —fDd^/dq shown in fig|| Although the vector field s(q) itself is 
isotropic its components are highly anisotropic fields. Again, let us consider the small 
angle approximation and t>i(q) oc si(q) component of the velocity vector field. The power 
spectrum P si (k) = kfP$(k) is anisotropic with considerably more power along k\ direction. 
Therefore the field si(q) has much smaller scale along q\ direction than along q<i direction. 
For an arbitrary direction hi = q%lq it is translated into considerable anisotropies between 
radial and transverse directions. A natural expectation suggests that the structure in 
redshift space is shifted to the segments of the zero contour v q (c[) = satisfying the 
convergence condition dv q /dq < marked by black dots on green contour in fig.||. On 
the inner side of these segments the radial component of the peculiar velocity is positive 
and therefore directed away from the center while on the outer side it is negative and 
thus directed toward the center resulting in shifting mass to these lines by the mapping to 
redshift space. 



o 
o 

N) 
O 
O 
-J 

o 
o 
o 



The progenitor of the structure in redshift space determined by eq. 3.2 is shown as 
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Figure 5: Lagrangian space. Contours fd qq (q) — 0.25,0.5,1 are shown in blue, green and red 
respectively. 



red regions in fig.|. The green contours show the boundaries of the regions <5(q) > 1 shown 
in fig.Q as green regions and also as green contours in figj2|. It is worth comparing the 
progenitors of the mappings to Eulerian (fig.§) and redshift space (fig.|7]). The circular 
pattern due to fDd gg term in eq. [T^ is quite obvious in fig.[7| but it is significantly less 
conspicuous than in fig.|| because the other two fields (Ii and I2) affecting the Jacobian are 
both statistically isotropic and homogeneous. The black dots mark the lines of convergence 
defined above and shown in fig|]. The convergence lines obviously correlate with the 
progenitor of the mapping to redshift space although the correlation is not perfect. 

The result of the mapping of the nonlinear progenitor (eq. |3.2| ) to redshift space is 
shown as red regions in fig|8]. If the linear progenitor was mapped to redshift space it would 
sit perfectly on top of the nonlinear structure, however it would remain significantly less 
connected than the nonlinear structure as the mapping is continuous. The figure illustrates 
a substantial change in the geometry caused by the mapping itself. The sizes of red filaments 
and voids between them are noticeably greater in redshift space than in Lagrangian space. 
Some smaller regions seen in Lagrangian space have merged into larger structures when 
mapped to redshift space. The difference between the structures in Eulerian (fig.§) and 
redshift space (fig.||) is quite substantial. 

Figure || also shows the convergence segments of the zero line of the radial velocity 
field demonstrated in fig|| Similarly to fig]?] the black dots are noticeably more numerous 
in red regions than between them. In this case the correlation is not entirely perfect either 
but it is quite obvious despite the fact that the black dots are plotted in Lagrangian space 
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Figure 6: Lagrangian space. Contours of fv q (q) > are shown in red, fv q (q) = in green and 
fv q (q) < in blue. The segments of convergence i.e., the parts of the contour v q (q) = satisfying 
the condition dv q /dq < are marked by black dots. 

while the red structure is plotted in redshift space. As we have seen in fig]!] the mapping 
to Eulerian space does not shift the progenitor much and the following mapping to redshift 
space does not move the points with v q = 0. In a few cases the separate segments of the 
convergence set seem to be aligned forming the filaments with total length of over 500 
h _1 Mpc. An example of this is the system of five segments that begins at distance about 
200 h. —1 Mpc toward lh from the center and continues clockwise down to distance ~150 
h _1 Mpc toward 4h from the center. 

The correlation between the structure in redshift space and the segments of convergence 
allows to use the lengths of the segment of convergence as predictors of the length of 
filaments. Figure ^ shows the segments of convergence i.e., the parts of line i> 3 (q) = 
where additional conditions dv q /dq q < is satisfied. These lines are attractors of the 
nonlinear structure in redshift space. Inspecting the figure one can easily identify the lines 
extended over 500 h -1 Mpc that suggests the existence of the filaments of similar sizes. The 
field Sq(q) is closely related to the initial potential since s q (q) = —d^/dq and therefore 
may be affected by an insufficient size of the box. Thus, the actual size of largest walls 
can be even greater than suggested by fig.|| in real universe. The full three-dimensional 
N-body simulations checking this prediction are currently on the way. 

3.3 Scales of Gaussian random fields 

In cosmology Gaussian fields are often characterized by several scales derived from moments 
< k 2 i >= (2-k 2 )- 1 f k 2 iP(k)k 2 dkk 2j of the power spectrum P(k) (e.g., @). The ratios of 
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Figure 7: Lagrangian space. The nonlinear progenitor of the structure in redshift space defined 



by the condition eq. ^2 is shown in red. Green contours show the regions defined by the linear 
condition <5(q) > 1. Black dots mark the segments of convergence. 

these moments generate a set of scales R 2 = 3 < k 2 i > / < k 2 ^ +l ^ >. Two particular scales 
R\ = R* (using notation of ref. ^jj ) and Rq are used the most often as the characteristics 
of linear density fields. The former is associated with the spatial distribution of peaks and 
is discussed in detail in p3| . The latter is referred to as the mean distance between crossing 
points of zero level |3(J. introduced the integral length scale L = £ _1 (0) J °° £(r)dr = 
(7r/2) Jq 00 kP(k)dk/ Jq 00 k 2 P(k)dk, the physical meaning of which is not quite obvious. An 
example of using other scales in cosmology is given in (3^]. The suggested model deals 
with the components of the deformation tensor djj(q) and displacement Sj(q) that are 
combined into radial components d qq (q) and s q (q). Here we consider the small angle 
approximation that fully described by d\\ and s\ fields. The scales can be characterized 
by the distribution of distances between successive zeros of a field along two axes (see 
e.g., pq| ). Figure |l(] shows the cumulative probability functions of zero crossing distances 
for four fields (5 = Ii,dn, v\ and <£) along two axes (qi and q2)- The density contrast 
field is isotropic and both cpfs are practically identical. The potential field also must be 
statistically isotropic but it depends strongly on long wave part of the power spectrum. In 
fact the variance cr\ formally logarithmically diverges at k — > 0. However, this is not a 



serious physical problem. Heavens |33| pointed out that in a fixed sample volume (e.g., the 
observable universe) the potential is well behaved, however we cannot determine the global 
zero level of the potential. He showed that by introducing A$(x) = 3>(x) — where is 
the value of the potential at the observer location, the problems related to the divergence 
of <t| are reduced to a weak dependence of the results on unknown value of ^o- Such fields 
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Figure 8: Rcdshift space. The large-scale structure is shown in red. The black lines are the parts 
of the contour v q (q) = w i (q)n i (q) = selected by the condition dv q /dq < 0. Two blue arcs mark 
the central 50 h^Mpc region where the adopted approximation (Z?|s|/<7 <C 1 ) is not very accurate. 

are known as locally homogeneous and locally isotopic fields and are discussed in Sec. 13 of 
ref. |fill |. This dependence of the potential on the long wave part of the spectrum explains 
huge variations in the tail of the cumulative function seen in fig.[K] (black solid lines). 
The scales along q2 are much greater than along q\ for both anisotropic fields d\\ and v\. 
The level lines of all fields are quite wiggly therefore there are quite a few short distances 
between crossing points. The overall sizes are probably better characterized by the tails of 
the distributions. Here are the median and the value marking the top 10% distances along 
qi and qi axes for four fields in h _1 Mpc: median (5(11, 11), d\\ (10, 18), s\ (21, 31) and $ 
(102, 84), top 10% 5 (30, 30), d n (21, 54), Sl (66, 117) and $ (397, 437). 
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4. Implications for three-dimensional universe 

The progenitor of the mapping to Eulerian space is determined by the interplay of three 
fields: one Gaussian 5 = I± and two non-Gaussian I2 and I3. It is easy to show that they do 
not correlate by the direct calculation of the corresponding integrals over the joint pdf (see 
Appendix). However the fact that two fields are non-Gaussian does not allow to declare 
them statistically independent and it makes the further analysis more challenging. 

Many features of two-dimensional model although not all are expected to be very 
similar in three dimensions. Nonlinear term proportional to I2 helps to fill the gaps between 
peaks of linear field 5 = ii(q) > 1 in 2D and therefore boosts the length of filaments in 
Lagrangian space up to 70 -100 h^Mpc. There is no reason to doubt that it works in 
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Figure 9: Lagrangian space (1024 h Mpc box). The black dots mark the segments of convergence 
i.e., the segments of contour f g (q) = where additional condition dv q /dq q < is satisfied. These 
lines are a kind of attractor of the structure in redshift space and their lengths are similar to the 
lengths of walls. 



a similar way in 3D. In addition, the term proportional I3 in eq. |2.3| is likely to boost 
the lengths of filaments even more. However, the further study is obviously needed both 
analytical and numerical. 

The mapping to Eulerian space itself does not change significantly the lengths of the 
filaments and does not displace them much from the initial positions in 2D. The mapping 
in 3D is expected to be similar in this respect. 

However, the mapping to redshift space differs in 3D significantly in a few respects. All 
three fields d qq , M qq (see eq. |2.9| ) and s q (see eq. 2J3) have scales in the radial direction much 
shorter than in two transverse directions. This causes the progenitor and the structure in 
redshift space to acquire a circular pattern in 2D while it gives rise to the origin of shell 
like walls in 3D. The walls in redshift space expected to be neither homogeneous no thin 
but their transverse sizes must be considerably greater than the radial thickness. Similarly 
to the two-dimensional case the greatness of 'great walls' comes from anisotropy of two 
Gaussian fields d qq and s q and in addition one non-Gaussian field, M qq derived from initial 
Gaussian field. As fig. [h] indicates the transverse coherence of both du and especially 
s\ field can easily exceed 100 h _1 Mpc. Two isotropic non-Gaussian fields I2 and ^3 also 
contribute to the formation of the structure in redshift space since the corresponding terms 
are boosted by factor 1 + / in the Jacobian (eq. |2.9| ). 

The convergence lines defined in 2D become the convergence surfaces in 3D, they are 
expected to provide a reasonable estimate of the transverse sizes of great walls. 



o 
o 

N) 
O 
O 
-J 

o 
o 
o 



- 15 - 



0.01 



5 

cjujqi) 
dn(q 2 ) 
vi(qi) 
Vi(q 2 ) 
4> 




10 



100 



1000 



Mpc/h 



Figure 10: Cumulative probability function of the zero-crossing distances along two axes gi and 
Q2- The top doted line marks the median values of the probability functions. 



5. Summary and discussion 

We propose a new theoretical model that describes the formation of the large-scale structure 
in redshift space. It naturally explains the sizes of the observed great walls and predicts that 
even greater walls will be discovered in future deeper and larger redshift surveys. Although 



based on the Zel'dovich approximation [26| the suggested model uses it in a new and very 
different way. Instead of traditional approaches that focused on the evaluation of density 
or positions and velocities of the particles we introduce the notion of the progenitor of the 
large-scale structure. It is defined as a set of regions in Lagrangian space that after mapping 
to Eulerian or redshift space (depending on the kind of mapping we study) end up in the 
regions where the particle trajectories have crossed and the mapping is no longer unique. 
Neglecting the thermal velocities and assuming smooth initial condition the boundaries 
of these regions in Eulerian or redshift space are caustics. In reality due to hierarchical 
clustering process the caustics do not form but the density of both matter and galaxies 
is still highest in these regions (eg |^|]). Even in the case of smooth initial conditions 
the evaluation of the density within the shell-crossing regions is not easy and we do not 
attempt to estimate it in this paper. Instead we utilize the Zel'dovich approximation for 
identification of the progenitor of the structure in Lagrangian space and then for mapping 
it to Eulerian or redshift space. In particular we demonstrate the role of various fields 
some of which are non-Gaussian or/and anisotropic. This is the most important factor in 
explanation of the lengths of filaments in Eulerian space. Non-Gaussian ^(q) an d ^(q) 
fields fill the gaps between peeks of Gaussian 5(q) = h(q) field significantly boosting the 
length of the progenitor of the filaments in Lagrangian space. It is worth stressing that 
the mapping of the progenitor to Eulerian space does not increase its lengths noticeably 



although it significantly affects its appearance. 

The mapping to redshift space has both the progenitor in Lagrangian space and the 
structure in redshift space significantly different than the mapping to Eulerian space. In 
both cases the thickness of the filaments in Eulerian space and walls in redshift space is 
probably the least accurate in this approximation. The thickness of filaments in Eulerian 
space is probably exaggerated [27]. The thickness of the walls in redshift space may be 
affected by nonlinear effects that change the velocities in the shell-crossing regions. How- 
ever, the estimate of the extent of the walls in the transverse direction is expected to be 
quite robust. 

Although both the suggested model and one suggested by || (BKP) are based on 
the Zel'dovich approximation (eq. there are a few radical differences between them. 
The current model can explain the sizes of both the structures in Eulerian and redshift 
space emphasizing that the structure in redshift space must have significantly greater sizes 
in transverse directions than in the line of sight direction and therefore the wall kind 
structures must be more conspicuous in redshift space. The BKP model discusses only 
the structure in Eulerian space. The BKP model uses the smoothing scale i?b defined 
by an inequality as{Rb) ^ 1 such that it corresponds to "the conceptual picture of a 
smooth stately large-scale single-stream flow" [|| while the current model uses the scale 
of nonlinearity and therefore assumes the crossing of trajectories and presence of multi- 
stream flows. Although it is not excluded that the further comparison with cosmological 
N-body simulations may slightly change the scale for better fit to N-body simulations it 
will not change the basic construction used by the model - the progenitor of the structure. 
The model defines the progenitor of the structure as a set of regions in Lagrangian space 
where the Jacobian of the mapping becomes negative indicating that these regions end 
up in multi-stream flow regions after mapping to Eulerian or redshift space. There is 
also one more principal difference. The BKP model uses only linear i.e., Gaussian fields 
and conditional correlations between them. This model utilizes both non-Gaussian fields 
(I2, 13, M qq ) one of which is anisotropic as well as anisotropic Gaussian fields d qq ,s q along 
with linear density contrast field. The displacement fields Si and s q also play important 
roles in the appearances of the structure in Eulerian and redshift spaces. The suggested 
approach allows easily to isolate the role of every field and thus helps to better understand 
the complex nonlinear dynamics of the structure formation. 

Further studies including the more detailed studies of the statistical properties of the 
fields involved in the mappings as well as a direct comparison with cosmological N-body 
simulations are obviously needed. Some work has already started and the results will be 
reported in the following papers. 
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A. Appendix: PDF of invariants /i,/ 2 and J 3 



From eq. 12 in |2£| it is easy to derive the joint distribution of three invariants if one notes 
that dl\ dl2 dl 3 = (Ai — A 2 )(Ai — As)(A 2 — As)d\i d\2 dX 3 (Aj and I\ are measured in units 
of as, I2 in units of <7 2 and ^3 in units of <r|) 

p(h,I 2 ,I 3 ) = -^exp(-3lf)F 2 (I 1 ,I 2 )F 3 (I 1 ,I 2 ,I 3 ), (A.l) 



where 



F 2 (h,h) = exp ^"^2^ , if h < if/3 and otherwise, 

F 3 (I X ,I 2 ,I 3 ) = 1, if 4 2) <I 3 < 4 1} and otherwise 

The constraints on ^3 reflect the fact that all three Aj are real which requires that A(I±, I 2 , 13) < 
0, where 

A(Ii, I 2 , h) = Till - 2/ 1 (9/ 2 - 2/2)73 + 7 2 2 (47 2 - if). (A.2) 
Therefore, 73 must lie between the roots of the equation A(7i, I2, 13) = which are 



r(l,2) 



1 r 



h 27 



/!(9/ 2 - 27?) ± 2(7? - 37 2 ) 3 / 2 l . (A.3) 



The constraint on 7 2 is a trivial consequence of relations between the eigen values Ai,A 2 
and A3 and the invariants: I\ = Ai + A 2 + A3, 7 2 = AiA 2 + A1A3 + A 2 A3 and 73 = AiA 2 A3. 
Thus, the distribution of 73 for given 7i and 7 2 is uniform in the range 1^ < 1% < I^\ 
The distribution of 7 2 for given I\ and 73 is exponential in the range — 00 < 7 2 < 7?/3 and 
the distribution of I± for given 7 2 and 73 is Gaussian. Integration over 73 gives the joint 
distribution of I± and 7 2 

P(h,h) = ^ exp (-37?) exp (^7 2 ) (7? - 37 2 ) 3 / 2 , (A.4) 

and then the integration over 7 2 gives the Gaussin distribution of the linear density 
contrast since 5 = DI\ in the linear regime. It is not difficult to check that the invariants 
do not correlate at one point: i.e.,< 7i7 2 >= 0, < /1/3 >= and < 7 2 73 >= 0. 
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